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ABSTRACT 

High-energy collisions of heavy ions provide a means to study QCD in a regime 
of high parton density, and may provide insight into its phase structure. Results 
from the four experiments at RHIC (BRAHMS, PHENIX, PHOBOS and STAR) 
are presented, and placed in context with the lower energy data from the AGS and 
SPS accelerators. The focus is on the insights these measurements provide into the 
time history of the collision process. Taken together, the data point to the creation 
of a deconfined state of matter that forms quickly, expands rapidly and freezes out 
suddenly. 
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1 Introduction 



The goal of high energy heavy ion physics is to study QCD in a regime of high 
temperature, high density, and large reaction volumes. The hope is to find conclusive 
evidence that QCD undergoes a phase transition at a critical temperature from 
a confined state, where quarks and gluons are bound in colorless hadron states, 
to a deconfined quark-gluon plasma (QGP), where quarks and gluons can explore 
volumes larger than the typical hadron radius (R ~ 1 fm). Lattice calculations, 
under a variety of assumptions (e.g. number of quark flavors or m s ), make the 
robust prediction that the degrees of freedom available to the system rise rapidly as 
it is heated through the critical temperature T c ~ 150 — 200 MeV JT[|, as shown in 
Fig. 0. This corresponds to an energy density of order e c ~ 1 — 2 GeV/fm 3 . There 
is great interest within the theoretical community whether these lattice predictions 
will be confirmed in experiments and, more generally, whether high temperature 
QCD can be used to make quantitative predictions. 
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Figure 1: Lattice predictions for the en- Figure 2: Theoretical expectations for the 
ergy density as a function ofT/T c . QCD phase diagram. 

While the theoretical questions are compelling in themselves, the exper- 
imental challenges of measuring the properties of heavy ion collisions in the labo- 
ratory imply a related set of questions that should be addressed before definitive 
theoretical conclusions can be made. The overarching problem is whether it is pos- 
sible to discern the various stages in the time evolution of the system, which evolves 
from the initial collisions to the final state hadrons in times on the order of 10 's of 
fm/c. To discover a QGP formed in the early stages of the collision, there must be 
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observables which maintain information about its properties throughout the expan- 
sion and cooling of the system and the final hadronization stage which forms many 
of the particles that register in the experimental apparatus. 

While measurables which address the early times directly (e.g. lepton and 
photon observables) are becoming available, the emphasis here will be on hadronic 
observables. The claim in this proceedings is that the evidence collected so far 
points to a picture of the collision dynamics which is consistent with the formation 
of a deconfined state that forms quickly, expands rapidly, and freezes out suddenly. 
There are also interesting comparisons with charged particle production in e + e~ 
collisions that provides some connections with perturbative QCD phenomenology. 

2 Colliders & Experiments 

High energy beams of heavy nuclei have been available at the Brookhaven AGS 
(Au+Au at the CMS energies per NN collision from ^s^n = 2.5 — 4.3 GeV). and 
the CERN SPS (Pb+Pb collisions from = 8-17.3 GeV). Since then, a new 

collider, RHIC, has begun colliding nuclei at higher energies (up to ^s^n = 200 
GeV) with substantial luminosities (of up to 10 26 sec _1 cm 2 ). The results shown here 
focus primarily on RHIC data from Au+Au collisions over a range of energies (^snn 
= 19.6, 56, 130 and 200 GeV). Most of the results are from the 130 GeV run, where 
the available luminosity limited the pt reach of many observables. However, early 
results from the 200 GeV data are shown here, and some preliminary results shown 
at the recent Quark Matter conference are also mentioned 0. 

There are four RHIC experiments, which can be classified into two groups: 
the "large" experiments, PHENIX and STAR, which are large-volume and large- 
acceptance general-purpose detectors, and the "small" experiments, BRAHMS and 
PHOBOS, which are have more limited acceptance but cover aspects of the collisions 
not addressed by the other experiments. STAR has a large TPC with coverage for 
charged hadrons in \r)\ < 1 and particle identification via dE/dx up to px ~ 1 
GeV. PHENIX has more limited hadron coverage (|r/| < .35) but a higher data rate, 
extensive particle identification for hadrons, photons, and electrons, and forward 
muon arms. PHOBOS has a small-acceptance silicon spectrometer near y ~ 1 and 
An acceptance for multiplicity measurements. BRAHMS has two movable multi- 
particle spectrometers, one at midrapidity, and one covering rapidities out to y ~ 4, 
offering a look at identified particle spectra in regions not covered by any of the other 
RHIC experiments. The four RHIC experiments complement each other in that they 
study many different aspects of heavy ion collisions, but all have hadron coverage 
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near mid-rapidity, which allows detailed cross-checks between the experiments. 
3 Initial Conditions 
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Figure 3: Ratios of particles to an- 
tiparticles vs. rapidity, from the 
BRAHMS experiment. 



Figure 4: Charged particle density at 
midrapidity, for e + e~ , pp and Au+Au col- 
lisions. 



The initial conditions of heavy ion collisions are determined primarily by 
the beam energy and the nuclear geometry. The impact parameter between the 
nuclei controls the overlap area of the contracted nuclei and thus the number of 
nucleons which participate via inelastic processes (N part ) and the number of bi- 
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nary collisions (N coU oc N,p art ) The beam energy determines the initial state parton 
density, which is dominated by gluons at low-a;, the available range in rapidity 
(Ay = 2 ln(-y/s/m p )), and the total nucleon-nucleon cross section. The combina- 
tion of these factors contribute to the amount of energy deposited by the incoming 
baryons in each nucleus as they lose energy while penetrating the other nucleus - a 
phenomenon called "proton stopping" This is typically measured by the net rapidity 
distribution of protons minus anti-protons. Recent results from BRAHMS show 
that while the ratio of protons to anti-protons is about 0.75 at mid-rapidity, it falls 
rapidly by y = 3, consistent with a rapid increase in the net proton density, which 
is expected to peak at a rapidity 2.5 units from ?/beam[f||- These measurements will 
be critical for understanding the very early collision stage. 

All of the RHIC experiments perform measurement at midrapidity, where 
hard processes are expected to contribute substantially to particle production. The 
energy density achieved in the collisions can be estimated either by measuring the 
transverse energy and using the Bjorken formula e = (dEx / 'dy) /tcR 2 t , with R the 
nuclear radius and r c the "formation time" of the system. PHENIX has measured 
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the transverse energy density at midrapidity in 130 GeV Au+Au collisions and 
extracted an energy density e = 4.6 GeV/fm 3 |J assuming at typical value for 
t ~ 1 fm/c. This is already well above lattice expectations for a QGP phase 
transition. They have also found the ratio Et/N^ to be constant for all N c h. Thus, 
one can use the ratio of midrapidity charged particle production at 130 and 200 GeV 
(R = 1.14 ± .05), shown in Fig. | as the particle production per participant pair ||, 
to infer an energy density of about e = 5.3 GeV/fm 3 at the highest RHIC energies. 




Figure 5: Centrality dependence of 
dN ch /dr]/(N part /2) compared with two- 
component and saturation models. 



Figure 6: Pseudorapidity distributions 
of Au+Au, pp and e + e~ at y/s=200 
GeV. The parton saturation prediction 
of Ref. J^/ is also shown. 



The centrality dependence of charged particle production, dN c h/drj/(N part /2), 
offers a means to study the evolution between pp collisions and central Au+Au. Two 
types of predictions @ have proven effective in describing the midrapidity data || 
shown in Fig. 5. The "two-component" model postulates that the total particle 
production stems from a linear combination of soft processes that scale with the 
number of participants, and hard processes that scale with the number of binary 
collisions. Thus, dN ch /drj = n pp ((l — X(s))N part /2 + X(s)N co n), where n pp is the 
measured multiplicity in pp collisions, and X(s) is the energy- dependent fraction of 
hard processes in pp collisions, (approximately 0.1 at RHIC). It is also expected 
that the energy released in the form of low-x gluons should create a parton density 
so high that the partons below a "saturation scale" Q s (i.e. of larger transverse 
size) recombine in order not to exceed a maximum value of order l/a s (s) ||. In this 
picture, particle production is determined mainly by Q s , which itself depends on the 
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Figure 7: Total charged particle produc- 
tion for a variety of strongly-interacting 
systems. Comparison with a fit to the 
e + e~ data is shown. 
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transverse density of partons (i.e. Q 2 oc A 1 / 3 ), dN c h/dr]/ '{N part /2) = l/a s (Q 2 ). As 
shown in Fig. |5], both of these models offer an efficient description of the experimen- 
tal data from all of the RHIC experiments. Saturation models are also able to offer 
a reasonable quantitative description of the full rapidity distribution, an example of 
which by Kharzeev and Levin || is shown as a dotted line in Fig. ||. 

The success of the saturation models in capturing basic features of particle 
production has a fascinating consequence. While most expectations of the "forma- 
tion time" t used in energy density estimates are around 1 fm/c, saturation models 
predict the formation time of the initial gluon state to be r s ~ h/Q s ~ .2 fm/c, 
which gives e ~ 18 GeV/fm 3 when plugged into the Bjorken formula for the 130 
GeV data. While neither estimate (5 GeV/fm 3 or 18 GeV/fm 3 ) should be taken too 
seriously at this point, these estimates show that both naive and more sophisticated 
estimates arrive at values well above the lattice expectations for QGP formation. 



4 Universal features of particle production 

We can gain insight into the relationship between Au+Au collisions and elementary 
nucleon-nucleon interactions by comparisons with e + e~ annihilations to hadrons. 
It has long been noted that the charged multiplicities in pp and e + e~ collisions are 
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similar if the energy taken away by the "leading baryons" is subtracted from yfs and 
then compared to e + e~ collisions at this v/s^f, suggesting a universal fragmentation 
mechanism 0. With new data from RHIC and LEP2, PHOBOS finds that Au+Au 
(divided by the number of participant pairs (N part /2)) and e + e~ collisions have 
similar pseudorapidity densities over a large range (see Fig. ^|) and thus a similar 
total multiplicity (TT|. Comparisons of the total multiplicity per participant pair 



(N C h) I {Np ar t/2) over a large range of CMS energies is shown in Fig. [7] and compared 
to e + e~ and pp/pp data. It is easiest to compare these data sets by dividing all of 
the values vs. \fs by a function that describes the e + e~ data [12). It is seen that the 



Au+Au and e + e~ data converge for */s > 20 GeV, while the pp data also follows 
the same trend if the effective energy ^/s^s = y/s/2 is used. 

The total multiplicity produced per participant pair in Au+Au collisions 
has also been measured by PHOBOS. It is shown in Fig. |8] to be constant over 
the range of centralities measured (N part > 65) for all three energies fllfl . This is 



reminiscent of "wounded-nucleon" scaling |13] but in all three cases the multiplicity 



is about 30-40% higher than pp at comparable energy, and comparable to e + e~. 
This adds an interesting perspective to the use of minijet and saturation approaches 
in describing the multiplicity at r\ = 0, as discussed in the previous section. 

This similarity in bulk particle production between Au+Au and e + e~ data 
for \fs > 20 GeV suggests that Au+Au are more efficient than NN interactions 
in transferring the incoming energy into particle production. This might simply 
indicate that the multiple collisions suffered per participant substantially reduce the 
leading particle effect. The drop below 20 GeV may be explained by the presence 
of a large baryon density in the final state fifM. While it is obvious that Au+Au 



collisions do not start with the same initial scale as e + e~, which is reflected in the 
much harder p? distributions in e + e~ above 2 GeV, it appears that the available 
energy, rather than the initial hard scale, has the stronger effect on the bulk particle 
production. This may be another manifestation of the universality of low- a; physics 
in strongly-interacting systems at high energy (which is discussed in the HERA 
context in Ref. fl5) ). The physics behind this may need to be addressed before 
claims to particular behavior for particle yields in heavy ions can be made in other 
contexts. 

5 Dynamics of the Parton Cascade 

If the energy density is indeed as large, and forms as early, as the success of the satu- 
ration descriptions seem to imply, then we might expect the density of scatterings to 
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Figure 9: Data from STAR and PHO- 
BOS on elliptic flow vs. collision central- 
ity compared to hydrodynamic predictions. 
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Figure 10: STAR data on V2 vs. px 
for different particles species, also com- 
pared with hydro calculations. 



lead to an early thermalization and possibly the development of hydrodynamic be- 



havior as the system evolves ||16|| . This expectation is borne out by measurements of 
the event-by-event azimuthal distributions, which are found to have a slightly ellip- 
tical shape at all energies from AGS to RHIC. The maximum value of V2, the second 
Fourier component of the azimuthal distribution, dN/d(j) oc 1 + 2t> 2 x cos(2(0 — \I/_r)) 
with the angle of the event reaction plane, has been found to increase with in- 
creasing beam energy. This is contrary to what one would expect if the dynamics 
were purely hadronic, since late stage rescattering would wash out any initial state 
anisotropy. Thus, the non-zero V2 has been interpreted as stemming from hydro- 
dynamic pressure gradients in the initial state, implying partonic thermalization, 
which is an essential prerequisite to QGP formation. 

The results from STAR and PHOBOS [|17|, are shown in Fig. |9] and 



compared with predictions from a range of hydro models (shown by solid bars). 
These results suggest that the more central events at RHIC are in fact consistent with 
hydrodynamic evolution. The models used here typically start with energy densities 
around 20 GeV/fm 3 evolving from an early time r < 1 fm, which is consistent with 
values extracted from the saturation approach. 

The importance of v 2 as an experimental observable comes from the quan- 
titative comparisons with theoretical predictions for observables relative to the re- 
action plane, which is defined as the axis where dN/d<p is at a maximum. Mea- 
surements of i>2 as a function of the transverse momentum |I7| characterize the 
modulation of dN/d(f> for particles of a given momentum relative to the reaction 
plane. As seen in Fig. [10], STAR has measured V2(pr) for different particle species 
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and found good agreement with hydro calculations out to px ~ 1 GeV, where their 
particle identification stops. However, two notable disagreements with hydro exist. 
The linear rise of v 2 with px stops at px ~ 2 GeV and remains constant out to larger 
Pt ]HJ • Also, the pseudorapidity distribution of v 2 drops rapidly away from midra- 
pidity, suggesting that the "boost-invariant" hydrodynamic behavior may only be 
found very near y = [TH . 



6 Thermal Description of the final Hadronic State 

As the system expands and cools, it eventually hadronizes into the final state parti- 
cles measured in the detectors. As these are the particles identified by most of the 
experiments, a large amount of data exists on the properties of this stage, only a 
subset of which will be shown in this section. The most pressing issue is whether 
the final state shows collective effects, which would be a signal of equilibration, as 
opposed to particle production simply proceeding via the available phase space. 

The most basic property of the system, measured by all four RHIC ex- 
periments, is the gradual approach of all of the anti-particle/particle ratios towards 
unity |19| , as shown in Fig. |ll|. The fact that p/p is less than unity directly indicates 



the presence of some fraction the primary baryons at mid-rapidity, an incredible fact 
considering the large initial baryon momentum. The kaons, while approaching one, 
also do not reach it at the top RHIC energy. These two ratios have been found to be 
highly correlated, suggesting that the produced hadrons are highly sensitive to the 
net baryon density. Still, RHIC is closer than ever before to creating the conditions 
found in the very early universe. 

The relationship between these ratios, and many others, can be understood 



by means of "thermal models" |2(| , which assume the system to be in chemical and 
thermal equilibrium up to the time when the momentum transfers are low enough 
that the flavor composition of the system is frozen. If only the ratios of parti- 
cle yields are used, then the available parameters are simply T (the temperature), 
fiB (the baryon chemical potential) and 7 S , which parametrizes the deviation from 
full strangeness equilibration. A volume (V) is needed if absolute yields are to be 
characterized as well. 

The application of these models to data from central collisions of heavy 
nuclei, as a function of -y/s, finds the simple result (shown in Fig. §) that all of the 
data seem to lie on a single contour on the (T, //g) space. There are two interpreta- 
tions of this behavior ||20||. Cleymans et al. explain this as freezeout occurring when 



the energy per particle (E) / (N) ~ 1 GeV. Braun-Munzinger and Stachel suggest 
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Figure 11: PHOBOS ratios of Figure 12: Thermal freezeout parameters shown as 
anti-particle to particle vs. y 7 ? a function of ^Js. 

that it occurs at a constant total baryon density of ns = .12/fm 3 . Both of these 
interpretations are consistent with the existing data. 

It is interesting that thermal models have also been applied to pp, pp and 
e + e~ data and find the similar result of T ~ 170 MeV (with fig = 0) ||21|| . This 
suggests that the heavy ion data above CERN energies approaches similar freezeout 
state as that for elementary systems. This may be another indication of the onset 
of universal behavior as the role of the net baryon density becomes less significant. 
However, some important remain between heavy ions and the elementary systems. 

One difference between nuclei and the more elementary systems is evident 
in the details of the momentum spectra. In elementary collisions (pp, e + e~), the 
spectra at low itlt show "m^ scaling" where the individual particle yields described 
by an exponential exp(— m^/T) with same slope for all particle species. In heavy 
ions, the slope at low tut increases with the particle mass in a way which can be 
fit by the form T e g = T Q + m(j3) 2 , where ((3) is a collective "radial flow" velocity. 
The energy dependence of the fit parameters has been extracted by Kaneta and 
Xu f22| and is shown in Fig. [12]. They observe a thermal freezeout temperature 
(T t h = 140 MeV) which is consistently lower than the chemical freezeout extracted 
from the particle ratios (T c h = 170 MeV). They also find that (/3) increases rapidly 
wity yfs but appears to saturate above ^Jsnn ~ 10 GeV. These data strongly imply 
collective behavior among the final state hadrons. 

Another important difference is seen in the total strangeness yield, which 
is "enhanced" by a factor of two in heavy ion collisions, as shown by the energy 
dependence of the "Wroblewski" factor (A s = 2(ss)/((uu) + (dd))) in Fig. |13| . 
While strangeness enhancement has long been a canonical signal of QGP formation, 
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Figure 13: Energy and system depen- 
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there exists a simple interpretation of this effect in terms of the volume over which 
strangeness can equilibrate. In other words, the system produced in nuclear collisions 



no longer seems to be forced to obey local strangeness conservation [23 ] . However, 
a recent analysis by NA49 [2J]|, characterizing the K + /tt + ratios as a function of 
centrality at CERN SPS energies, has shown that the effect may not be related 
precisely to the reaction volume. Instead, this ratio scales with the fraction of 
participants that have collided more than once, shifting the responsibility for what 
appears to be equilibration in the final state to some property of the initial collision 
geometry. 

Both chemical and thermal freezeout is addressed in the the model by 
Broniowski and Florkowski [Q . They assume that chemical and thermal freezeout 
occur at the same time. By fitting to the available 7T, K and p spectra from PHENIX 
and STAR, and incorporating transverse expansion, this model is able to describe 
all measured particle species, and was even successful in predicting particles like the 
K* and Q. The success in understanding the K* yields is very interesting since it has 
a lifetime comparable to the lifetime of the system itself (cr ~ 4 fm). This suggests 
that there is little, if any, subsequent interactions of the K* decay products, lending 
credibility to the assumption of the single freezeout time. This model is similar in 
concept to so-called "blast wave" fits to experimental data, which also assumes that 
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hadronization occurs on a thin space-time shell. The success of these two approaches 
suggests that the hadronization process is quite sudden. 

7 Hard Probes of the Early Stages 

Hard processes occur at early times (At ~ h/AE) and thus offer a means to probe 
the early stages of the collision. In principle, jet cross sections can still be cal- 
culated within a pQCD framework as they are for hadron-hadron interactions, by 
considering the structure of the nucleons in the colliding nuclei and assuming they 
interact via the standard pQCD matrix elements. Naively, one would expect that 
the rate of hard processes should scale with the number of binary collisions, since 
the timescale of hard processes is so short that incoming nucleons will be sensitive 
to all of the nucleons in its path. However, there is large body of theoretical work 
suggesting that high-energy partons should lose energy in a deconfmed medium via 
color brehmsstrahlung - a phenomenon called "jet quenching". Jets traversing a 
hadronic medium should suffer no induced radiation and thus be unmodified. 

At RHIC energies, there are two impediments to directly measuring jets. 
First of all, the typical energy is not high enough to create well-collimated cones of 
particles. Secondly, a central event always creates a large background of soft particles 
which will obscure a jet signal within a typical jet cone. Thus, measurements of 
quenching focus primarily on the high-p-^ part of the inclusive hadron spectrum, to 
look for modifications to the measured jet fragmentation functions. 

First results on particle spectra at high px in heavy ion collisions were 
shown by the PHENIX and STAR collaborations only a few months after the end 
of the 130 GeV RHIC running in 2000 |26|| . Both collaborations presented the mod- 
ification of the inclusive hadron spectra (both of all charged particles and identified 
7r 0, s) in central Au+Au collisions by a ratio comparing it with spectra from pp colli- 
sions measured by the UA1 collaboration, interpolated to 130 GeV and then divided 
by the number of binary collisions. 

dN/dp T {Au + Au) 
AA ~ N coll x dN/dp T (pp) { } 

The PHENIX results of Raa vs. Pt for charged hadrons and neutral pi- 
ons in 130 GeV Au+Au collisions are shown in Fig. [DJ At low-pr, the Raa is 
1/6, as expected for wounded nucleon scaling (since N co u/(N part /2) ~ 6 in central 
Au+Au). However, while Raa increases with p T , as might be expected for the 
gradual dominance of hard processes, at 2 GeV the rise stops, with Raa saturat- 
ing or decreasing, depending on the particle type, well below unity. STAR results 
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on Raa show a marked decrease above 2 GeV. Preliminary results at 200 GeV 
from both experiments confirm this decrease and find that above 5-6 GeV, Raa be- 
comes approximately constant at ~ 1/6 out to 9 GeV (PHENIX all charged and tc°) 
and ~ 1/3-1/2 out to 11 GeV (STAR all charged). These results are in qualitative 
agreement with the jet quenching hypothesis and theoretical descriptions are rapidly 
becoming available. 

While large energy loss effects were expected for the nearly-massless light 
quarks, heavy quarks are expected to radiate far less. This is because of the "dead- 
cone" effect, where radiation is suppressed at angles less than 9 < m q /E, where m q 
is the mass of the heavy quark, and E is its energy [ 27j . This expectation has been 
tested by PHENIX using measurements of single electrons and positrons in central 
130 GeV Au+Au collisions. They observe a significant excess of electrons above 
the expected hadronic and photonic backgrounds and attribute it to the presence 
of open charm, thus extracting an open charm cross section per binary collision 
dec = 380/i6 ± 200(sys) ± 60(stat) p8fl . This is consistent with an extrapolation of 
FNAL and ISR measurements using a PYTHIA calculation tuned to reproduce the 



lower energy results, as shown in Fig. 16 



Thus, there seems to be no violation of collision scaling in the charm sector, 
consistent with no energy loss, as expected from the dead-cone effect. However, it 
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must be kept in mind that this measurement has large systematic errors, which 
are being carefully addressed in the higher-statistics 200 GeV data. It should also 
be noted that the electrons have not been directly tagged as coming from charm 
decays, e.g. by measuring the characteristic decay length of the charm mesons. 
Future upgrades to the RHIC experiments seek to address this issue definitively, 
which is a crucial piece of information both for understanding jet quenching and as 
a baseline for J/ty suppression. 



8 Ultra Peripheral Interactions 
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Figure 17: STAR results on coherent p° pro- 
duction in 130 GeV ultraperipheral Au+Au 
collisions. 



The colliding gold beams at 
RHIC have also opened up opportu- 
nities of studying the physics of co- 
herent meson production. The highly 
compressed ions act as sources of very 
strong Coulomb fields as well as pomerons. 
Thus, at large impact parameter, one 
expects photon- photon or photon- 
pomeron interactions, as well as vir- 
tual photon exchange. Vector mesons 
can be produced coherently by a vir- 
tual photon, emitted by one nucleus, fluctuating into a quark-antiquark pair and 
scattering elastically off of the other nucleus, provided they satisfy the conditions 
p T < 7th/ Ra (90 MeV) and p\\ < h?ij/Ra (6 GeV) where 7 is the lorentz boost of 
the nucleus (70 at RHIC 130 GeV). The p° production cross section can be predicted 
using Glauber extrapolations of p° photoproduction data to be 350 mb [f29] . 

Using a combination of special triggers and minimum-bias samples, STAR 

has extracted a strong coherent p° peak both in collisions where the gold nuclei are 

intact, and where one or more neutrons is dissociated from each nucleus via coulomb 

excitation (shown in Fig. [H]). They extract an exclusive p° cross section of 410 ± 

190 ± 100 mb PP| . The agreement between these measurements and the Glauber 

predictions suggests that p° production and Coulomb excitation are independent 

processes, and thus factorize. These measurements thus point to new opportunities 

in the physics of strong fields. 
9 Summary of Results 



With the new RHIC data, systematic data now exists for heavy ion collisions as 
a function of \/s over several orders of magnitude and as a function of impact 
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parameter. These data test the interplay between hard and soft processes in a large- 
volume system where nucleons are struck multiple times. The data is consistent with 
creating a deconfined state (jet quenching) that forms quickly (saturation models), 
expands rapidly (radial and elliptic flow) and freezes out suddenly (single freezeout 
and blast wave fits). There are also intriguing connections with particle production 
in elementary systems, which point to the role of the energy available for particle 
production on the features of the final state. Many in this field are optimistic that 
the careful understanding of this experimental data may lead to the theoretical 
breakthroughs that will connect these complex systems to the fundamental lattice 
predictions. 
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